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ABSTRACT 

We study the mid-infrared properties of 1315 spectroscopically confirmed members in eight massive 
(M V i r > 5 x 10 14 M Q ) galaxy clusters covering the redshift range from 0.02 to 0.83. The selected clusters 
all have deep Spitzer MIPS 24/xm observations, Hubble and ground-based photometry, and extensive 
redshift catalogs. We observe for the first time an increase in the fraction of cluster galaxies with mid- 
infrared star formation rates higher than 5 M© yr _1 from 3% at z = 0.02 to 13% at z = 0.83 (Rp ^ 
IMpc). This increase is reproduced even when considering only the most massive members (M* ^ 
4 x 10 10 M Q ). The 24/im observations reveal stronger evolution in the fraction of blue/star- forming 
cluster galaxies than color-selected samples: the number of dusty, strongly star- forming cluster galaxies 
increases with redshift, and combining these with the optically-defined Butcher-Oemler members 
[A(B — V) < —0.2] doubles the total fraction of blue/star-forming galaxies in the inner Mpc of the 
clusters to ~ 23% at z = 0.83. These results, the first of our Spitzer/MIPS Infra-Red Cluster Survey 
(SMIRCS), support earlier studies indicating the increase in star-forming members is driven by cluster 
assembly and galaxy infall, as is expected in the framework of hierarchical formation. 
Subject headings: galaxies: clusters: general - galaxies: evolution - galaxies: fundamental parameters 



1. INTRODUCTION 

iButcher fe Oemlerl (|1978l 11984ft observed that galaxy 
clusters at intermediate redshift have a higher fraction 
of members with blue optical colors than clusters in 
the local universe, thus providing a key piece of ev- 
idence supporting galaxy evolution. This increase in 
blue members with redshift, named the Butcher-Oemler 
(BO) effect, was intensely debated for two decade s 
(e.g. iMathieu fc Spinradl fl98lt iDressler fc Guml fl98l) . 
However, multiple optical studies based on spectro- 
scopic observations have since confirmed the increase 
in blue, star- forming galaxies in higher redshift clus- 
ters (e.g.lCouch fc SharSelll987tlCaldwell fc Rosi 19971: 
iFisher et alJl!998t lEllingson et al.ll2001h . and found that 
BO galaxies reveal signs of recent and ongoing star for- 
mation. The paramount question now is have we seen 
only the tip of the iceberg? 

Most studies of star-forming galaxies in clus- 
ters rel y on rest- frame ultraviolet or optical trac- 
ers (e.g. iBalogh et all 119981: iPoggianti et al.ll2006f ). but 
UV/optical tracers can suffer severely from dust obscu- 
ration, especially when star formation is concentrated 
in the nuclear regions (|KennicuttJ Il998f h For exam- 
ple, ultraluminous infrared galaxies have SF rates of 
> 1000M Q , yet many ULIRGs fail to even be detecte d 
at UV and optical wavelengths (e.g. iHouck et al.l f2005). 
Although corrections for dust attenuation are possible, 
reliable estimates of SF rates cannot be achi eved solely 
using rest-frame U V/optical observations (jBelll 120021 : 
iCardiel et alll200lh . 

A substantially more robust method of determining to- 
tal SF rates is with mid-infrared (MIR) imaging. The 
first MIR imaging of galaxy clusters at intermediate 
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redshifts was tak en with ISO's ISOCAM camera, and 
IDuc et all (j2002ft found that at least 90% of the star 
formation was hidden at optical wavelengths. The first 
handful of galaxy clusters observed with the MIPS cam- 
era on the Spitzer Space Telescope (SST) have also re- 
vealed strong dust-obscured star formation (jGeach et al.l 
2006; iMarcillac et alll2007t iBai et al]l200l . 

It remains unclear as to what causes the increase in 
star-forming galaxy cluster members. Detailed morpho- 
logical studies of blue galaxies [defined as having A(B — 
V) < -0.2] 3 with the Hubble Space Telescope (HST) 
find that most are disk systems similar to those in lo - 
cal clusters (e.g. IDressler et alJll994t ICouch et al.lll994ft : 
past studies also fi nd that many show s i gns of inter- 
actions or mergers (lLaver~ fc Henry 1988; L averv et al.l 
fl99l lCoucLlt^fl99i lOemler et al1ll997n . More re- 
cently, studies indicate that galaxy infall is a viable ex- 
planation for the significant numbers of blue galaxies 
and their dist urbed morphologies in intermediate redshift 
clust e rs (e.g. Ivan Dokkum et al.l fl~998bt lEllingson et alJ 
120011: iTran et alJ 12005), a scenario suppor ted by hierar- 
chical clustering models (jKauffmannl 1995). In this case, 
galaxy clusters that are accreting a significant number 
of new members should have a higher fraction of star- 
forming galaxies, especially at higher redshifts when the 
amount of activity was enhanced also in the field. 

Here we present the first comprehensive study of 
SST/MIPS 24/itm imaging of galaxies that are spec- 
troscopically confirmed members of eight massive 
(M V ir >5x 1O 14 M ) X-ray luminous clusters spanning 
a wide redshift range (0.02 < z < 0.83). After presenting 
the data in §2, we focus our analysis in §3 and §4 on 
the evolution of star- forming members with redshift. A 
cosmology with (iJ ,^A/,^A) = (70 km s -1 , 0.3, 0.7) 
is assumed throughout the paper; at z = 0.83, the 

3 A(B — V) is the color offset from the red sequence fit to the 
cluster ellipticals. 



2 



Saintonge, Tran & Holden 



look-back time is ~ 7 Gyr. 



2. DATA 

We have assembled a data set of eight galaxy clusters 
at 0.02 < z < 0.83 that have a total of 1315 spectroscopi- 
cally confirmed members. The core of our sample is com- 
posed of five clusters spanning the entire redshift range 
with large spectroscopic membership, uniform multi- 
filter optical photometry and deep SST/MIPS imaging 4 . 
For the part of the analysis that does not depend on 
rest-frame (B — V) color, we fold into the sample three 
additional clusters: A bell 1689 for whi ch MIR data from 
ISOCAM is available (jDuc et al.ll2002l ). and CL0024 and 
MS0451, both of w hich have extensive redshift catalogs 
(jMoran et all [2005) and MIPS observations. Observa- 
tional details for all clusters are in Table [TJ 

2.1. Optical Photometry and Spectroscopy 

Th e optical photometr y for the five main clusters is 
from iHolden et alj (|2007l hereafter H07) where magni- 
tudes and colors were derived from Sersic models fitted to 
HST/WFPC2 images (MS1358, MS2053, and RXJ0152), 
HST/ACS images (MS1054), and SDSS mosaics for 
Coma. The conversion to rest-frame values is done 
by in terpolating between the passbands (|Blakeslee et al.l 

2006) and has errors of ~ 0.02 mag. The mass-to-light 
ratios (M/Lb) and stellar masses were calculated using 
the relation between rest-frame (B — V) color and M/Lb; 
see H07 for details and a discussion on the associated er- 
rors. 

2.2. MIPS 24fim Imaging 

All MIPS data sets were retrieved from the Spitzer 
public archive. Individual frames were corrected with 
scan mirror position-depend ent flats and then mosai ced 
with the MOPEX software (jMakovoz fc Khanll2005T l to 
a pixel size of 1.2" 5 . Integration times (ti nt ) and 
background levels (F& g ) in these mosaics are given 
in Tabled) Photometry was performed with APEX 
(jMakovoz fc Marleaul [2005) using a 3"-diameter aper- 
ture, and an aperture correction of 9.49 as given in the 
MIPS data handbook. A small aperture is necessary 
to avoid contamination in the deep and crowded clus- 
ter fields. The fluxes are consistent with results from 
PSF-fitting photometry with scatter from a 1:1 relation 
in the range of 15-25 /iJy. 

To estimate the completeness of each MIPS catalog, 
we added to the mosaics artificial sources modeled on 
the PSF. To avoid overcrowding, we simulated 30 signals 
at once, and repeated the process 30 times for each clus- 
ter (the 50% completeness limits, F 50 %, are presented in 
Table [TJ. Finally, the MIPS sources were mat ched with 
the o ptical catalogs using a 2" search radius (|Bai et alJ 

2007) . From randomization of the MIPS coordinates, we 
estimate the rate of false identification to be 7 ±4%, and 
little dependency of this error rate on redshift or color is 
observed. 

4 This work is based on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. 

5 The instrumental pixels are 2.55"in size, but the finer sampling 
helps in improving the characterization of the PSF. 



2.2.1. Star formation rates 

Star formation rates are based on the 24/im fluxes. 
First, the total infrared luminosity (-Fs-iooo/jm) of each 
galaxy was determined using a fa mily of infrared spectra l 
energy distributions (SEDs) from lDale fc Heloul (pTjol ). 
We choose the range of SEDs that are representative of 
the g alaxies in the Spi tzer Infrared Nearby Galaxies Sur- 
vey ([Dale et al.l l2007h , and at each redshift adopt the 
median conversion factor from i^^m to Fs-iooo/xm given 
by these models. At 0.4 < z < 0.6, the error due to 
the adopted conversion factor is ~ 20%, but the error 
increases to a factor of 1.5-2.0 at lower and higher red- 
shifts. For the parts of our analysis that are sensitive 
to the SF rates, we take these errors into account. As 
a check, we note that our total infrar ed lumino s ities i n 
MS1054 agree well with the values in iBai et all (|2007h . 
The conversion from total infrar ed luminos i ties t o star 
formation rates is done following iKennicu tt (1998). 

We assume that the emission at 24^m is due to star 
formation but it could also be due to dust-enshrouded 
active galactic nuclei (AGN). However, in comparing the 
X-ray and 24/im detections, only one cluster galaxy (in 
RXJ0152) is detected in both and rejected. While the 
AGN fraction in clust ers seems to increase with redshift 
(|Eastman et al.l |2007[). the estimated AG N fraction is 
only 2% at z ~ 0.6. Ijohnson et all ([20031 ) also find evi- 
dence that at z ~ 0.8, any excess X-ray AGN are located 
at R > 1 Mpc whereas we focus on the central Mpc of 
each cluster. Although we cannot completely rule out 
possible contamination by weak obscured AGN, we have 
excluded X-ray AGN and thus assume that the galaxies 
det ected by MIPS a re pow ered by dusty star formation; 
see iMarcillac et al.1 (|2007l ) for a more detailed argument 
on why this is a reasonable assumption. 

3. RESULTS 

3.1. Color-Magnitude diagrams 

Figure [TJ presents the color-magnitude diagrams of the 
five main clusters with photometry from H07. Because 
the MIPS sensitivity varies from cluster to cluster, we 
apply a SF rate limit of 5 M Q yr _1 . The first immediate 
observation is that the number of strongly star-forming 
galaxies increases significantly with redshift. Using a 
field galaxy sample drawn from the same photometric 
and spectroscopic catalogs, we estimate a possible field 
contamination at z = 0.83 to be ~ 8% (i.e. no more than 
one galaxy per cluster). In Figure [TJ the dotted lines rep- 
resent the original color criterion for BO galaxies. The 
ratio of the number of cluster galaxies with MIR SF rate 
> 5M yr _1 above this color cut to the number of blue 
galaxies (A(B — V) < —0.2) increases with redshift. 

3.2. The Mid-Infrared Butcher- Oemler effect 

For each cluster, we compute and plot in Figure [2] the 
fraction of confirmed star-forming cluster members af- 
ter selecting by rest-frame B-band magnitude (Mb < 
— 19.5), cluster-centric distance 6 , and MIR star forma- 
tion rate (^ 5Mq yr _1 ). The errors on fsF,MiPS repre- 
sent the range that can be produced by taking the min- 
imum and maximum conversion factors from f24/jm to 

6 While the optical observations generally extend to R > 1.5 
Mpc, the MIPS imaging for MS1358 only extends to ~ 1 Mpc 
(~ 50 - 60% of r200 for these clusters). 
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TABLE 1 
Properties of selected clusters 



Name 


Coords [J2000] 


z range a 


Lx h 






^24° 


tint 


F bg 


-^50% 








(10 44 erg s- 1 ) 








(s pix ) 


(MJy sr" 1 ) 


(M© yr" 


Coma 


125935.7+275734 


0.013-0.033 


9.0 + 0.2 


244 


63 


134 (2) 


73 


32.8 


0.02 


Abell 1689 


131129.5-012017 


0.17-0.22 


21.4 + 1.0 


81 


52 


12 (2) 


f 




1.4* 


MS 1358+62 


135950.4+623103 


0.315-0.342 


10.2 ± 0.7 


171 


73 


21 (3) 


2700 s 


20.6 


0.75 


CL 0024+17 


002635.7+170943 


0.373-0.402 


2.9 + 0.1 


205 


51 


11 (6) 


2700 g 


48.5 


1.48 


MS 0451-03 


045410.9-030107 


0.52-0.56 


21.0 + 0.4 


242 


38 


8 (5) 


2700 g 


35.0 


3.34 


MS 2053-04 


205621.3-043751 


0.57-0.60 


6.5 + 0.4 


85 


13 


15 (8) 


1950 


35.2 


5.04 


MS 1054-03 


105700.0-033736 


0.80-0.86 


16.4 + 0.8 


112 


75 


13 (8) 


3600s 


47.4 


4.54 


RX J0152-13 


015243.9-135719 


0.81-0.87 


18.6 ± 1.9 


117 


61 


19 (8) 


3600S 


31.9 


3.05 



a Cluster members selected within this redshift range, as in H07 (Coma, MS1358, MS2053, MS1054, RXJ0152), IDuc et all J2002D 
(A1689) a nd IMoran et alJ 1120071) (CL0024 an d MS0451). b B o lomet ric X-ray l uminosities from H0 7 (Coma, MS1358, MS2053, MS1054, 
RXJ0152), IBardeau et all J2007TI (Abelll689l. IDonahue et al.l ( 119991) (MS0451), IZhang et all (120051) ( CL0024). C Total number of spectro- 
scopically confirmed members (magnitude-limited selection). Rcdshifts arc taken from Bcijcrsbcrgcn (2003) Due ct al. (2002) Fisher ct al. 
d 19981) [Moran et al.l d2005f ). IMoran et al.l (120071) iTran et al.l (120051) ITran et al.l d2007DlDe"marco et all ( I2005D . respectively." 1 Number of con- 
firmed members within 1 Mpc of the cluster center and brighter than Mb — —19.5 + 5 log h. c Number of MIPS det ections in the clu ster; 
() is the number of galaxies within N s with SF rates ^ 5 Mq yr~ . We are using ISOCAM mid-IR data from IDuc et al.l 112002!) for 
A1689. g Over the central 5'x5'of the MIPS image. 
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Fig. 1. — Rest-frame color-magnitude diagrams for spectroscopi- 
cally confirmed members in the main cluster sample. Filled red cir- 
cles are MIPS detections with SF rates 5 Mq yr — 1 (where all clus- 
ters are better than 50% complete). The larger symbols represent 
galaxies with log 1(1 (M*) 10.6. The rest-frame B-band magnitude 
has been corrected for passiv e luminosity evolution, as d etermined 
from the fundamental plane (Ivan Dokkum et al. 1998a). The ver- 
tical dashed line is the rest-frame B-band magnitude selection limit 
of -19.5. The solid diagonal li ne is the best fit to the red sequence 
galaxies, adopting the slope of van Dokkum ct al. (1998b), and the 
dotted line denotes A(B — V) = —0.2 mag; only galaxies below the 
dotted line would be part of standard BO sample. 

-Fs-iooo^m instead of a single average value for each clus- 
ter, and by varying the different selection thresholds by 
amounts comparable to the errors on each of these pa- 
rameters. 

Figure [5] shows that the fraction of galaxies in clusters 
with MIR SF rates ^ 5M Q yr _1 steadily climbs from 
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Fig. 2. — Fraction of confirmed cluster galaxies that are star- 
forming as revealed by the MIPS 24/^m observations. Are consid- 
ered only members with MIR SF rates ^ 5 Mq yr — 1 that are 
brighter than Mg = —19.5 and located within 1 Mpc of the clus- 
ter centers (filled circles) and 500 kpc (open squares). The points 
for the two z ~ 0.83 clusters are offset slightly in z for clarity. 

~ 3% locally to ~ 13% at z = 0.83. Because H07 showed 
that a cluster's morphological composition can vary de- 
pending on whether members are selected by mass or 
by luminosity, we apply an additional stellar mass cut 
of log 10 (M„) ^ 10.6 (Fig. [3]). The mass cut is only ap- 
plied to the five main clusters for which uniform photom- 
etry and thus stellar masses are available; the remaining 
three clusters are shown only as upper limits. While 
the mass cut attenuates the increase in fraction of star- 
forming members, it does not completely suppress the 
trend. Thus the MIR BO effect is not due to an increase 
in the fraction of faint, low- mass members temporarily 
brightened by strong star formation. 

4. DISCUSSION 

Having established an increase in the fraction of MIR- 
detected galaxies from z ~ to z ~ 0.8, we stress that 
optical studies are likely underestimating the increase 
in star- forming cluster galaxies with redshift. As seen 
in Fig. [TJ an increasing number of strong dust-obscured 
star-forming members appear on or near the red sequence 
at higher redshifts; these are not included in traditional 
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Fig. 3. — As in Fig. [2] the fraction of confirmed star-forming clus- 
ter galaxies (MIR SF rate> 5M yr _1 , M B < -19.5, R < 1 Mpc), 
but now with the additional stellar mass cut-off of log 1 Q(M«) ^ 
10.6 for the five main clusters (filled circles). Stellar masses are 
not available for the remaining three clusters; they are shown as 
upper limits (open circles). 

color-selected BO studies. The late- type morphologies 
of these members supports our intepretation of dusty 
star forma tion and red colo rs due to extinction (see also 
A901/902: IWorf et alJl2005h . 

Using the standard BO definition of A(B - V) < -0.2, 
the fraction of blue galaxies with Mp < —19.5 and 
Rp < 1 Mpc at z ~ 0.8 is ~ 11%; however, including 
the red, massive, star-forming members raises the total 
fraction of blue/star- forming members to ~ 23%. We 
note that for the five main clusters, the increase in the 
blue/star- forming fraction due to these red, star- forming 
members is {1.1, 1.2, 1.3, 1.7} at z ={0.02, 0.33, 0.59, 
0.83}, i.e. the relative importance of including red, dusty 
star- forming members increases with redshift. 

Is this increase linked to galaxy infall? In Figure [2j 
both CL0024 - 0.4) and MS2053 (z ~ 0.6) are above 
the general trend established by the other six clusters. 
Both CL0024 and MS2053 have enhanced star formation 
compared to other clusters at similar redshift, and both 
have bimodal redshift dist ributions. CL0024 is made of 
two colliding subclusters (Cz oske et al.l [2002) , and has 
an unusually larg e number o f luminous infrare d galaxies 
(|Coiaet al.ll2005h . Similarly, EaneFaD $2QQB) conclude 
from that MS2053 has a significant number (> 25%) of 
infalling galaxies; these members tend to be blue and 
star-forming. Both CL0024 and MS2053 are accreting 
a large number of new members and have high frac- 
tions of dusty star-forming galaxies. We speculate that 



the increase in star-forming members reflects the recent 
accretion of new members, i.e. galaxy infall, and that 
such events are more frequent at higher redshift due to 
the process of cluster assembly ([Ellingson et al.l 120011 : 
iTran et al.ll2005t lLoh et all 120081 ). As further evidence 
of this, 80% of the MlPS-detected galaxies in the z ~ 0.8 
clusters are more than 700 kpc from the cluster cores in 
projected distance and thus the MIR Butcher-Oemler ef- 
fect is significantly altered by only considering the inner 
500 kpc of the clusters (open symbols in Figf5]). 



5. SUMMARY 

We present the first comprehensive study of 
SST/MIPS 24/im observations for seven massive, 
X-ray luminous galaxy clusters spanning a wide redshift 
range (0.02 < z < 0.83). Uniform photometry, high 
resolution HST imaging, and extensive redshift catalogs 
enable us to measure the fraction of members with 
strong, dust-obscured star formation. The fraction 
of cluster galaxies with MIR star formation rates 
^ 5M Q yr" 1 increases from 3% in Coma to ~ 13% in 
clusters at z = 0.83, and this trend is evident in both 
luminosity (Mp < —19.5) and mass-selected samples 
(M* > 4 x 1O 1O M ). 

Optically-based studies increasingly underestimate the 
total amount of star formation in cluster galaxies with 
redshift because many of these dusty red star-forming 
members are missed in color-selected samples. These 
tend to be late-type galaxies that are red because of dust 
extinction which disguises their high levels of obscured 
star formation (> 5M Q yr _1 ). Defining the SF frac- 
tion to include both optically blue and red, but MIPS- 
detected members doubles the fraction at z = 0.83 from 
~ 11% to ~ 23% (R P < IMpc). 

Lastly, our study indicates that the BO effect and the 
increase in obscured star-forming members are linked to 
galaxy infall: 80% of the MIR-detected members at z ~ 
0.8 are outside the cluster cores (Rp > 0.7Mpc), and the 
two clusters at z < 0.8 that are accreting a substantial 
number of new members also have an enhanced fraction 
of galaxies with MIR SF rates > 5M Q yr" 1 . 
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